Self-Organized Criticality Below The Glass Transition 
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We obtain evidence that the dynamics of glassy systems below the glass transition is characterized 
by self-organized criticality. Using molecular dynamics simulations of a model glass-former we 
identify clusters of cooperatively jumping particles. We find string-like clusters whose size is power- 
law distributed not only close to T c but for all temperatures below T c , indicating self-organized 
criticality which we interpret as a freezing in of critical behavior. 
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Although the relaxation dynamics of glass-forming liq- 
uids has been studied for decades, there arc still many 
unresolved questions [1]. Especially for the dynamics 
of the glass out of equilibrium it is still an open and 
hotly debated question what characterizes the relaxation. 
We present in this Letter work on the dynamics below 
the glass transition where we focus on cooperative mo- 
tion. Cooperative rearranging regions have been stud- 
ied mostly above the glass transition and are the basis 
of Adam-Gibbs theory [2]. Above the glass transition 
two kinds of cooperative motion have been identified: (i) 
string-like motion [3-5] where of the order of ten parti- 
cles move along a Conga-line and where each particle is 
significantly more mobile than an average particle, and 
(ii) very cooperative motion where of the order of 40 par- 
ticles participate and where each particle undergoes only 
a small displacement [6]. 

To study cooperative motion below the glass transition 
we use molecular dynamics simulations, which have the 
advantage of providing us with the microscopic informa- 
tion of every particle's position at all times. Using these 
particles trajectories we first search each simulation run 
for jump events where a particle jumps out of its cage 
of neighbors. Then wc identify clusters of cooperatively 
jumping particles, i.e. jump events which are correlated 
in space and time. We find that the cluster size distri- 
bution follows a power law independent of details of the 
cluster definition. Furthermore, we find string-like clus- 
ters as they have been found above the glass transition. 

A similar cluster definition for the system under study 
but at a temperature slightly above the glass transition 
has also revealed a power law distribution [3]. Such dis- 
tributions are a signature of criticality, such as the clus- 
ter size distribution in percolation theory at the critical 
point [7]. However, contrary to these simulations and 
percolation theory, we find a power law not only close 
to a critical point but for all temperatures below T c that 
we have investigated (see Fig. 1). We thus find a type of 
self-organized criticality. 
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FIG. 1: Distribution P(s) of cluster sizes of simultaneously 
jumping particles for temperatures T = 0.15 — 0.43. The in- 
set shows P(s) at T = 0.43 where the jump events have been 
divided into five equal time windows. The black solid line repre- 
sents the first window. 



Our data are consistent with the following scenario: a 
glass cooled down to T c develops critical behavior, and 
then upon further cooling the criticality remains frozen 
in. Mode-coupling theory for glasses predicts the devel- 
opment of critical fluctuations [8] when T c is approached 
from above and a recent extension of mode-coupling the- 
ory to temperatures below the transition has predicted 
such a freezing in of critical behavior [9]. The signa- 
ture of criticality in these theories is the power law of 
a two-time correlation function [9-11]. Our data, while 
consistent with the mode-coupling freezing in scenario, 
find criticality in a spatial structure rather than a relax- 
ation exponent. We believe this represents the first direct 
observation in glasses of self-organized criticality, that is, 
criticality for all temperatures below T c . 

Our system is a well-studied binary Lennard- Jones 
(LJ) mixture of 800 A and 200 B particles. We refer 
the reader for details of the model to [12] and for de- 



2 



tails of the molecular dynamics simulations to [13]. Pre- 
vious simulations have shown that this system exhibits 
the main features of glass-forming liquids and is thus a 
good simple model for glass-formers [12]. For present-day 
computer simulations the system falls out of equilibrium 
in the vicinity of T g = 0.435 (in reduced LJ units) [12]. 
Whereas Donati et al. [3] studied cooperative motion of 
this system above the glass transition, we study here the 
same system but below the glass transition at 10 temper- 
atures ranging from 0.15 to 0.43. We use 10 independent, 
well-equilibrated configurations at T = 0.5 and then in- 
stantly quench the system to the desired temperature, 
e.g. T = 0.15. After an (NVT) run of 2000 time units we 
then run the (NVE) production run for 10 5 time units. 

For the definition of jump events we use the trajec- 
tory r n (t) of each particle n and take time averages 
over 800 time units to obtain its thermal fluctuation a n 
and average positions r^(ti) at times ti — 800(7 — 0.5) 
where I = 1,2,..., 125. We define a particle n to un- 
dergo a jump if its change in average position Af^ = 
\r^(ti) - r^(t;_ 4 )| satisfies Ar^ > \/20ct„ [14]. We thus 
identify for the whole simulation run all jump events 
{n,li,{r n )i,lf,(r n )f} of jumping particles n, jumping 
from average position (r„)j at time t^, the time asso- 
ciated with bin k, to average position (r„)/ at time ti f 
[15]. 

To address the question of cooperative motion we in- 
vestigate how these single particle jump events are corre- 
lated in time and space. To identify correlations in time, 
we group the jump events according to the bin index Zj. 
We thus obtain Ni simultaneously jumping particles for 
each time bin I. To investigate how these ./V; particles are 
spatially correlated, we identify clusters where particles 
n and m arc defined to be neighbors (and therefore mem- 
bers of the same cluster) if their distance |(r„)j — (r m )i\ 
is smaller than the position of the first minimum r m ; n 
of the corresponding radial pair distribution function of 
the complete system (r m ; n = 1.4 for AA, 1.2 for AB and 
1.07 for BB independent of temperature). This analysis 
gives us for each time Ki > distinct clusters. The clus- 
ters are numbered by k — 1, 2, . . . , Ki and we denote by 
Wj ; fc the set of particle labels composing the kth cluster 

K, 

in time bin I with Ni t k particles (i.e. ^l,k — A;). We 

fe=i 

now look at the size distribution Pis) of all clusters Afi t k, 
with s being the number of cluster members, i.e. 

^) = £X>(*,aWE^ (i) 

I fe=l 'I 

where 5(x,y) is the Kronecker delta function. 

In Fig. 1 we show a log-log plot of P(s) for tempera- 
tures T = 0.15 — 0.43. We observe essentially a straight 
line, i.e. a power law P(s) ~ s~ r , indicating scale invari- 
ance. The exponent for different temperatures ranges 
from t = 2.2 — 2.7 with the trend of r increasing with 
temperature. The average is r = 2.50 ± 0.05. A power 
law of P(s) has also been found in simulations of the same 



binary Lennard-Jones system slightly above but close to 
T c with r = 1.86 [3]. In percolation theory the size distri- 
bution n s (where P(s) = n s j J2 S n s) follows also a power 
law where the mean-field exponent r = 2.5 and in three 
dimensions r = 2.2 [7]. However, in these simulations 
and in percolation theory the power law occurs only at 
T c . In contrast, we find a power law for all temperatures 
(see Fig. 1), i.e. we find scaling invariance for the whole 
temperature range below T c . As mentioned above, this 
is consistent with the scenario of critical behavior being 
frozen in, and remaining for all temperatures below T c . 

This scale-invariance independent of a control param- 
eter is usually found in systems with self-organized crit- 
ically [16] and is, to our knowledge, a new phenomenon 
for the cluster size distribution of structural glass form- 
ers. The occurrence of self-organized criticality is usu- 
ally associated with being out of equilibrium and having 
widely separated time scales. Our system is consistent 
with these requirements: single particle jumps take of 
the order of 10 time units, the time between successive 
jumps is of the order of 30000 [13], and the equilibrium 
relaxation time is significantly longer than the simulation 
run [12]. 

Due to the importance of implications for relaxation 
dynamics below the glass transition, the question arises 
if our results for the cluster size distribution are specific 
to details of the analysis. We therefore investigate how 
robust P(s) is with respect to different definitions. Since 
the system is out of equilibrium, we first ask the ques- 
tion if the jump dynamics changes during the simulation 
run. We therefore divide the whole simulation run into 
five time windows of equal size and determine P(s) with 
Eq.(l) by restricting the sum over I accordingly. The 
inset of Fig. 1 shows the resulting P(s) for T = 0.43. 
We find that in the first time window the most collec- 
tive processes occur: up to approximately 250 particles 
jump simultaneously and spatially correlated. The dis- 
tribution P(s) however seems to follow a power law with 
t = 2.5 ±0.1 not only independent of temperature but 
also independent of the time window (i.e. waiting time). 
We find the same waiting time independence also for all 
other temperatures. 

To further check the sensitivity on details of the power 
law of P(s) we next modify the definition of a cluster. 
Whereas usually cluster connections are defined in space 
(and therefore avalanche-like correlations are tested in 
space), we generalize now the definition of a cluster by 
treating space and time similarly (and thus allowing for 
avalanche-like correlations in time also). Instead of re- 
quiring as before that two jump events a and (3 occur 
simultaneously (If = if) we define extended clusters by 
allowing two jump events to occur at neighboring time 
bins (i.e. |AZ| < 1). We show in Fig. 2 results for 
two different definitions of extended clusters (I & II). 
The difference between these definitions is due to the 
usage of time and position before or after the jump. 
We define two jump events {n a , If, (r n )f, IJ , (r„) J} 

and \mP , if , (r m )f , If, (r m )^} to be connected if 
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def.I: \lf - if | < 1 and |(r n )f - (r m )f | < r min 
dcf.II: (|Z? - Zj| < 1 and |(r n )f - (r TO )J| < r mln ) or 
- /f | < 1 and \(r n )J - (r m )f | < r ml „) 
As before for simultaneously jumping particles, defini- 
tions I and II result in a cluster size distribution which 
follows a power law, however, with a smaller exponent of 
t = 2.13 ± 0.04 and r = 2.22 ± 0.04 for clusters I and II 
respectively. As shown in Fig. 2 for T = 0.30 and 0.43 
we find again a power law for all temperatures (similar 
results are obtained for other temperatures). 
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FIG. 2: Distribution of cluster size for extended clusters I & II 
for temperatures T = 0.30 and 0.43. 



dence of (z(s)) and therefore an additional average over 
simulation runs at different temperatures has been in- 
cluded in Fig. 3. The comparison with an ideal string and 
the most compact cluster (sphere) indicates that both the 
clusters of simultaneously jumping particles as well as the 
extended clusters are string-like. This is similar to the 
results of cooperative motion above the glass transition 
[3-5] and below the glass transition [17]. 

Our results are consistent with the following scenario: 

above the critical temperature T c string-like clusters are 
i i i i I i i i i I i i i i I i i i r 
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FIG. 3: Average coordination number z within a cluster as a 
function of cluster size s. Error bars are of the order of A(z) = 
0.5. 



To further investigate these clusters we next character- 
ize their geometric shape directly via coordination num- 
bers (instead of displacement-displacement correlations 
as in [3]). We determine for each cluster A/j,fe the average 
coordination number 



(2) 



where z n is the number of neighboring particles m € A/"/, ft 
of particle n. Fig. 3 shows the average 



K, 



= E E s ( s > /EE ^ s < w.*) ( 3 ) 
i fc=i ' i fe=i 

as a function of s. We observe no temperature depen- 



found. Close to T c the distribution of cluster sizes fol- 
lows a power law. Below the glass transition this critical 
behavior gets frozen in. Independent of details of the 
cluster definition and independent of waiting time, we 
find string-like clusters with a cluster size distribution 
which follows a power law for all investigated tempera- 
tures. We expect this self-organized criticality to occur 
also for other glasses out of equilibrium and leave the test 
thereof for future work. 
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